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ABSTRACT 



We perform an analysis of two-body non leptonic decays of B and B s mesons 
in the factorization approximation. We make use of the semileptonic decay 
amplitudes previously calculated on the basis of an effective lagrangian sat- 
isfying chiral and heavy quark symmetries and including spin 1 resonances. 
Exclusive semileptonic -D-decay data are used as experimental input. Our re- 
sults compare favorably with data, whenever they are available and indicate 
a positive value for the ratio of non leptonic coefficients 02/01, similarly to 
previous studies. 



Chiral and heavy quark symmetries [JTJ] |2j have been used to construct an effective 
lagrangian for light and heavy mesons including their effective weak interactions ||. This 
approach has been extended to include the nonet of l - vector meson resonances, 
coupled through the so-called hidden gauge symmetry, so that one is able to describe 
the interactions among heavy Qq a mesons (Q heavy quark, a = 1,2,3 for u,d,s), the 
pseudoscalar and the light vector meson resonances. In reference [[5], we performed an 
extensive analysis of heavy meson semileptonic decays into final states containing light 
mesons. The main result of the analysis was that, using experimental information on 
D — ► 7ilb>i, D — ► Klui and D K*lv\ as inputs, together with the hypothesis of a simple 
pole behaviour of the form factors in the whole region of q 2 , one can make predictions 
on the form factors describing the semileptonic decays of the B and B s mesons into light 
mesons. 

In the present letter we wish to apply these results to the evaluation of two body non 
leptonic decays || of B and B s within the factorization approximation. Before going 
into the details of our work, let us discuss the meaning and the possible limits of the 
approximation. 



1 Non leptonic amplitudes 

Let us first write down the weak non leptonic hamiltonian we shall employ in our calcu- 
lations. To be specific let us consider AB = 1, AS = AC = transitions (other relevant 
cases can be handled similarly). We take J7| 

r nl = % £ W [Ci(M + cmoI] (i.i) 

V " q=u,c 

where the four local operators Of, 0%, 0\ and 0\ are 

0\ = ^(1-75)^7m(1-75K (1.2) 
O u 2 = ^(1-75)^7m(1-7 5 K (1-3) 



and 0\ and 0\ are obtained by the change u — > c. In \L.2\ and |1.3| i and j are colour 



indices, whereas in |TTT| Ci(/x) and C^/i) are QCD coefficients computed at the scale 
\i ps m b . We have neglected, in writing eq. [O], the contribution of the QCD penguin 
operators O3, . . . Oq since their Wilson coefficients C3, . . . Cg are numerically very small as 
compared to C% and C2 || ; while these operators could be relevant in some processes where 
the current-current operators 0\ and O2 are not active or strongly Cabibbo suppressed, 
in the applications we shall consider here their contribution would never exceed 10% of 
the result making their neglect reasonable. 

The values of the coefficients C\ and C2 have been computed in || beyond the leading- 
log approximation; for mj = 4.8 GeV, m t = 150 GeV and A A ^ S = 250 MeV one finds 
C x = 1.133 and C 2 = -0.291. 

In order to discuss the factorization approximation let us be specific and consider, 
e.g., the decay B° — > 7r + 7r~. To compute < 7y + ty^\Hnl\B > we have to consider the two 
hadronic matrix elements < 7r + 7r _ |Oi|S° > and < 7t + tt~\02\B° >■ As for the first one, 
the factorization procedure amounts to write (T^ = 7^(1 — 75)): 

< TT+vr-idifi >=< 7r + |uiy>|/5° >< 7riJr^|o > (1.4) 
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This approximation can be given some theoretical justification on the basis of 1/N C ex- 



pansion || or color transparency |10[ . Moreover, from a phenomenological point of 
view there are indications that in the realm of B physics, factorization roughly works 
satisfactorily JT1J. As to the second operator, we can write 



< 7i + 7i-\O 2 \B >=< it + it-\^-0 1 + 2 \B° > (1.5) 



with 

2 = uT»T a bdT (1.6) 
In |1.6| T a (a = 1, ... 8) are colour matrices, normalized according to TrT a T b = ^S ab 



and the sum over repeated indices is understood. In the naive factorization, inserting 
the vacuum, one finds that the 2 operator, which contains coloured currents, does not 
contribute and therefore the contribution of H^l to the amplitude will be given only by 
the 0\ operator with a multiplicative coefficient a\ — C±+ ^-C 2 . However, as discussed in 
and [ 13 1 , the rule of discarding the operators with coloured currents while applying the 



vacuum saturation is ambiguous and unjustified. This is already one reason, among many 
others, to make the choice ]12| fl4| to treat a\ and the analogous parameter a 2 , multiplying 



2 in the naive factorization approximation, a 2 = C 2 + C\/N C} as free parameters. This 
will be our attitude in this letter as well. 

Let us recall that the analysis of D non leptonic decays leads to the empirical finding 
cii m C\,a 2 « C 2 | I2| . There has been some recent theoretical effort to understand 



the empirical rule of "discarding the 1/N C term ", which has shown that, at least in 
some channels (e.g. B° — > D + n~ ) and within a certain kinematical approximation (i.e. 
Mb,Md — > oo while their difference is kept fixed) there is a tendency to cancellation 
between the two terms appearing in |1.5| JT5] ]. Since these analyses are far from having 
reached a definite conclusion, we are justified in our choice, especially because from our 
analysis, as we shall discuss below, a positive value of a 2 appears favourite, a result that 
seems general jnj and would indicate that the positive 1/N C term cannot be neglected. 

In such a type of calculation one has to be aware of the many uncertainties and 
incompleteness heavily bearing down at practically all the levels of the procedure. Such 
uncertainties and incompleteness are unfortunately in a way or in another common to all 
calculations of such complex phenomenon as non leptonic decays, and like everybody else 
we cannot avoid them. 

The standard approach to non leptonic decays consists, as we have already described, 
in incorporating short distance effects into an effective hamiltonian constructed through 
operator product expansion and use of renormalization group, while leaving long distance 
QCD effects within the hadronic matrix elements. Already the construction of the effective 
hamiltonian contains uncertainties, particularly in cases where many different scales are 
present, such as for transitions from b to c. One has to add to this the general impossibility 
of obtaining a scale independent result when the matrix elements of the scale-dependent 
effective hamiltonian are taken within some hadronic model with usually unrelated scale 
dependence. One can only hope that, within a class of processes, one may choose a 
suitable scale at which this procedure approximately applies. 

The idea behind the factorization approximation is that hadronization appears only 
after the amplitude takes the form of a product of matrix elements of quark currents which 
are singlets in color, thus allowing for approximate deductions from semileptonic processes. 
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Different kinematical situations may suggest that factorization may apply better to some 
non leptonic processes rather than to others. For instance one intuitively expects that it 
may work better when a color singlet current directly produces an energetic meson easily 
escaping interaction with the other quarks. Independent of this, various other effects such 
as more or less strong role of long-distance contributions, including final state interactions 
effects, of small annihilation terms, more or less sensitivity to choice of the scale, etc., 
may suggest that the uniform simultaneous application of the factorization approximation 
to different processes must be subject to detailed qualifications and may at the end turn 
out to be incorrect. Unfortunately at the present stage of the subject one is forced to 
first collect informations by comparing an admittedly very rough procedure to the data 
available. 

We also note that the effective lagrangian approach can be used to write down the 



non leptonic hamiltonian in terms of the heavy and light meson fields however in 
this way one can make prediction only for a small region of the phase space where the 
light mesons are soft. 

According to what already said, in the factorization approximation, two body non 
leptonic decays of B and B s mesons are obtained by the semileptonic matrix elements of 
the weak currents between different mesons. A suitable form is 

< P(p')\V»\B(p) >= [ip + pT + Ml ~ 2 Ml q^iq 2 ) - M| ~ M V F (g 2 ) (1-7) 



<V(e,p') 



V - A")\B(j>) > 
i(M B + My) 



2V(q 2 



Mb + My 
e* ■ q 



e e uPaPj3 



q 



+ ie ■ q 



(M B + My) 
2M V 



q, 

(p + P')n 



Mq 2 



Ml -My 



A2{q 2 ) 



q z 



q^A (q 



where P is a light pseudoscalar meson and V a light vector meson, q = p — p' 

My 



Ao(0) 



2Mi 



2Mi 



(I- 



(1.9) 



v 



and Fi(0) = F (0). 



For the q 2 dependence of all the form factors we have assumed a simple pole formula 
F(q 2 ) = F(0)/(1 — q 2 /Mp) with the pole mass Mp given by the lowest lying meson with 
the appropriate quantum numbers (J p = + for F , 1~ for F\ and V, 1 + for A 1 and A 2 , 
0- for A ). 

The values of the form factors at q 2 = 0, as given by the study of the semileptonic 
decays as performed in 0, are reported in Table I. The errors in the Table follow from 
the experimental semileptonic exclusive D decays, that we used as inputs. In Table II we 
report the values of the pole masses we employ. Let us explicitly stress that the values 
reported in Table I were obtained in the leading order of the 1/Mq expansion of the 
heavy chiral current. In || we reported also the results of a fit obtained introducing 
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mass corrections to the leptonic decay constant ratio /b//d = \J M D /M B (modulo logs) 
predicted by the heavy quark effective theory. However, as discussed below (see also [Q). 
a few non leptonic B decays seem to disagree with this latter solution and therefore we 
choose to work everywhere at the leading order in 1/Mq expansion. 

As for the current matrix elements between a pseudoscalar meson P or vector meson 
V we take: 

< 0\A»\P(p) > = ifppf 
<0\V\V(p,e)> = f v M v e^ (1.10) 

In Table III we present the values of fp and fy used in computing the rates. We observe 
that recent CLEO data on the decay D s — > /iu point to a rather large value of fo a '■ 



fu a = 315 ± 46 MeV [17]. On the other hand all theoretical approaches indicate for fjj 
and fo s smaller values; for example QCD sum rules [[0| give f& = 223 ± 27MeV and 
fr) s = 277 ± 13MeV; a recent QCD sum rules analysis for the ratio /d 4 //d gives the 



value ~ 1.1 [20] and similar results are obtained from Lattice QCD or potential models 
(for a review see |2"I|]). In view of these results we have assumed for f D and f Ds the 
values reported in Table III that are compatible with both CLEO data and theoretical 
expectations. 



2 Numerical results 

Let us begin our numerical analysis by discussing a class of decays that depend only on 
the parameter 02- Recent data from CLEO Collaboration Jl7j allow for a determination 
of this parameter. From BR(B -> KJ/ijj) = (0.10 ± 0.016) x 10" 2 , BR(B -> K* J/ij>) = 
(0.19 ± 0.036) x 10~ 2 and BR(B -> Kip(2s)) = (0.10 ± 0.036) x 10" 2 we obtain 

|a 2 |~0.27 (2.1) 

Using this value we can compute the results of Table IVa that contain two body non 
leptonic B° decays depending only on a 2 and on current matrix elements between B and 
light mesons. We have used the recent determination for V u b arising from the CLEO 
analysis of the end-point of lepton spectrum in semileptonic inclusive B decays |22 ]: 
Vub/Vcb = 0.075. We take V cb = 0.040. 

Let us now study the coefficient a\. The best way to determine it is to consider B° 
decays into D* + tt~, D* + p~, D + tt~, D + p~ final states. In order to use the experimental 
data we need an input for current matrix elements between B and D, D* states. 

In [[J we did not consider B — >• D and B — ► D* transitions in the heavy quark effective 
theory; this subject has been recently actively investigated by several authors and we rely 
on their work to compute the corresponding non leptonic decay rates. The relevant matrix 
elements at leading order in 1/Mq are [|| 



< D(v')\V"(0)\B(v) >= yjM B M D £(w)[v + v'f (2.2) 
< D*(e,v')\{V» - A^ l )\B{v) >= 

^M B M D £(w) [-e^ xpT e* x v p v' T +i(l + v v')e^ - i(e* ■ v)v'»] (2.3) 



where w = v ■ v ' = (M§ + M|> — q 2 ) ^MbMd, v and v' are the meson velocities, and 



£(w) is the Isgur-Wise form factor. £(w) has been computed by QCD sum rules f2~3 



|24J , potential models, and estimated phenomenologically fll4j . We shall take for it the 
expression |25| 

f M = Tl — exp 
\l + w/ 



"(2p 2 



2.4) 



w + 1 

which reproduces rather well the semileptonic data[^6| with p ~ 1.19, V c b = 0.04, tb = 
1.4ps. We stress once again that we have chosen to work at the leading order in 1/Mq 
which is why we have not introduced the non leading form factors discussed e.g. in J2 - 



From the new CLEO data [|HJ BR(B° -> D+tt") = (2.2 ± 0.5) x 10~ 3 , BR(B° -> 
D*+ir-) = (2.7 ± 0.6) x lO- 3 , BR(B° -> L>+p~) = (6.2 ± 1.4) x 10~ 3 and BR(B° -> 
= (7.4 ± 1.8) x lO- 3 one gets 

|ai| ~ 1.0 (2.5) 

For completeness we have reported in Table IVb the results for non leptonic decays de- 
pending only on a>i and on the Isgur-Wise form factor £(w). 

In Table IVc we use the fitted value of ai together with the weak matrix elements 
computed in our model to compute additional two body decays depending on a±. We 
observe explicitly that the recent CLEO result [f£J BR(B° -> = (1.3±g;| ± 0.2) x 

10~ 5 is in agreement with our outcome. We stress that this agreement depends on our 
choice of the "scaling" solution for the semileptonic matrix elements, according to our 
previous discussion. The other solution discussed in ref. U appears to be excluded by the 
present data for this decay channel. 

In Table V we consider B~ decays. Let us observe that they depend on the relative 
sign between a% and a2 which has not been fixed yet. The new CLEO data [T7| BR(B~ 
D°n-) = (4.7±0.6) x 10" 3 , BR(B~ -> D*°tt-) = (5.0 ±1.0) x 10~ 3 , BR(B~ -> D°p~) = 
(10.7±1.9)xl0- 3 and BR(B~ -> D*°p-) = (14.1 ±3.7) x HT 3 , allow to conclude that the 
ratio a 2 /ai is positive. Clearly this result depends on the relative phase of the hadronic 
matrix elements. Our choice is the "natural" one, i.e. we assume (analogously to \T2\) 



that for a decay B — > M\M 2 {M\ and M 2 scalar or vector mesons) the phase is the one 
determined under the assumption of spin and flavour symmetry in the meson spectrum. 
Of course these symmetries are (even badly) broken in many decays; one should therefore 
be aware of the possibility to have a different phase between the two terms in such cases. 

Our results for B~ decays are reported in Table V, while in Table VI we compute 
analogous results for the B s non leptonic decay channels. In doing the calculations we 
have neglected the mixing rj-i]', and we have assumed ideal mixing between uj and <fi. 
Moreover we have not taken into account final state interactions, whose effect we cannot 



evaluate at the moment (in D decays, as shown in fT2" |, their effect was in some cases 



rather significant), and all effects of mixing and CP violation. 

We also note that our reported theoretical errors in Tables IV, V and VI only reflect 
the uncertainties in the weak couplings reported in Table I. Since the form factor Aq at 
q 2 = has an uncertainty of about 100% (see Table I), the BR's containing this coupling 
have a large error, as one can see in some entries of the Tables IV- VI and should therefore 
be taken with care. 
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Let us finally comment on the decays with two charmed mesons in the final state 
(see Tables IVb, Va and VIb). Our results are in agreement with data from both CLEO 



and ARGUS collaboration |30j |L7j. According to the discussion of ref. |pl| , however, 
the correct value of the a\ coefficient to be used in connection with the the factorization 
hypothesis and the Heavy Quark Effective Theory is |oi| ~ 1.45. Indeed, as discussed 
in PH| , since the light quarks do not carry large momenta in these decays, the running 
strong coupling constant should be computed at a scale lower than m c , which results 
in a value for the |ai| coefficient significantly larger than the one used in Tables IV- 
VI. Using experimental data on decays into £)(*)£)(*) [30J JX7j ; f or example (BR(B° — > 



D + D~) = (0.6 ± 0.45) x 10~ 2 ) one can see that the above mentioned large value of |ai| 
is incompatible with the data. It would be interesting to analyze if such disagreement 
reflects a breakdown of the factorization approximation or the presence of non leading 
effects in 1/m.Q approximation. 

In conclusion, we have performed, in the factorization approximation, an analysis 
of two body non leptonic decays of the B and B s mesons; our study has been based 
on semileptonic amplitudes obtained by an effective lagrangian having chiral and heavy 
quark symmetries and has employed semileptonic exclusive D decays as an input. Our 
results are in agreement with the experimental data, whenever they are available, and 
indicate, similarly to other analyses, a positive value for the ratio of the non leptonic 
coefficients a 2 /ai. Our results represent, in our opinion, a preliminary indication that B 
semileptonic decays to light mesons can be related to the analogous D decays without 
major violations of the heavy quark flavour symmetry. 
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Tables Captions 



Table I Form factors at zero momentum transfer for B — > P and B — > V semileptonic 
transitions (see ref.||). Indicated errors are from the experimental inputs. The form 
factors A\ and are relatively well known, while Aq has a relative error ~ 100%. 

Table II Pole masses for different states. Units are GeV. 

Table III Values of leptonic decay constants. For particles that can couple to different 
weak currents, the corresponding weak current is indicated within brackets. For <p 
and us we assume ideal mixing, while rj is the pure octet component. 

Table IV Predicted widths and branching ratios for B° decays. We use in the Tables 
t~b° = t b s = Tb- = 14 x 10~ 13 s, V u b = 0.003, V c b = 0.04. The quoted errors come 
only from the uncertainties of the form factors of Table I 

Table V Predicted widths and branching ratios for B~ decays. 

Table VI Predicted widths and branching ratios for B° s decays. 
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Table I 


Decay 


Fi — F 


V 


A, 


A 2 


A 3 = A 



B - 
B - 
B - 
B s 
B s 
B - 
B - 
B - 
B s 
B, 



K 
V 

* 7] 

> K 
P ± 

K* 

* K* 
,0 



0.49 ±0.12 
0.53 ±0.12 
0.49 ±0.12 
0.52 ±0.12 
0.52 ±0.12 



0.62 ±0.12 
0.62 ±0.12 
0.61 ±0.12 
0.64 ±0.12 
0.62 ±0.12 



0.21 ±0.02 
0.21 ±0.02 
0.20 ±0.02 
0.20 ±0.02 
0.19 ±0.02 



0.20 ±0.08 
0.20 ±0.08 
0.20 ±0.08 
0.21 ±0.08 
0.20 ±0.08 



0.24 ±0.24 
0.24 ±0.24 
0.20 ±0.20 
0.17±0.21 
0.17±0.18 



Table II 


State J F 


0" 


\- 


0+ 


1+ 


dc 


1.87 


2.01 


2.47 


2.42 


sc 


1.97 


2.11 


2.60 


2.53 


ub 


5.27 


5.32 


5.78 


5.71 


sb 


5.38 


5.43 


5.89 


5.82 



Table III 


Particle 


fp 


Particle 


fv 




132 


P ± 


221 


K 


162 


K* 


221 


D 


240 


D* 


240 


D s 


270 


D* s 


270 


rj(uu/dd) 


54 




409 


rj(ss) 


-108 





221 






UJ 


156 
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